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ABSTRACT 


Computer-aided maps have become a necessity in 
present day Military applications. Such a system, besi- 
des creating, storing and displaying the map, should 
have efficient query processing capabilities. In Mili- 
tary applications, answers to the graphical queries may 
be approximated with errors which are permissible within 
the specified limits. Curves, describing features on the 
map, are defined by sets of discrete digitized points. 
These curves can be approximated with different resolu- 
tions using Strip Tree concept. Based on this concept, 
a package has been developed. The package performs 
graphical query processing pertaining to the features on 
the map, which are represented by curves, using hierar- 
chical Strip Tree representation. A syntax-oriented 
command language has been designed for this purpose. 
Using an existing package for creating the base map, the 
package developed answers the queries, performing geome- 
trical computations on Strip Trees. 



CHAPTER I 


INTRODUCTION 


1.1 General 

With advancement of weapon technology , faster 
mobility and latest means of communications , modern war- 
fare has become highly complex and technical in nature. 
Guided missiles, armoured personnel carrier vehicles , 
computerised communication systems using satellite are 
some of the examples. It is very vital for a commander 
to have every possible information about strength and 
deployment of his forces, their fighting capability and 
morale, type of weapon systems, guns and fire power avai- 
lable, means of communication with his senior and subordi' 
nate commanders, administrative and replenishment facili- 
ties available, nature of ground and its topography where 
he is going to fight his battles etc. etc. 

1.1.1 Sources of Military Information 

As discussed in preceding paragraph, the nature of 
information required for military planning is very large 
and diverse in nature. It is therefore not possible to 
obtain it from one single source. 
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In an operational headquarters, such information 
is collected by different departments and stored in the 
form of reports, tables, charts, maps etc. For example, 
the department dealing with military intelligence 
collects the information from various sources and hands 
it over to the operation branch. Then this information 
may be used to update the tables or update the deploy- 
ment of forces shown on the map or sand models. Subse- 
quently when a commander has to plan an operation the 
relevant information is extracted from these different 
sources as and when required. 

1.2 Role of Maps and Sand Models 

Knowledge of ground, its topography and character- 
istics is one of the important feature for any military 
operation. The tactics for battle will depend greatly 
on this information. Some of the examples highlighting 
this aspect are: 

(i) Hill tops and high grounds are very .suitable 
location* for deploying forces in defensive 
role, 

(ii) An undulating ground like ravines will not be 
suitable for an operation involving tanks. 

(iii) Information regarding roads and tracks is vital 
to plan movement of troops and maintaining a 
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supply line. 

(iv) Rivers , Nalas and Canals are obstacles in the 
advance. The information about bridges, the 
fordability, current flow etc. are important 
for planning of crossing these obstacles by 
different means , 

Based upon this information, as well as other infor- 
mation of tactical nature, the commanders decide the 
deployment of their forces in detail. The positions of 
operational headquarters, which maintain the command and 
control of battle, various bodies of troops, guns, tanks 
etc. are decided. However, all these deployments are very 
dynamic in nature and undergo changes very frequently. 

The two commonly available methods to project the updated 
locations on the ground ares 

(i) Maps and Sketches 

(ii) Sand Models. 

1.2.1 Maps and Sketches 

Maps are generally the standard documents published 
by Survey of India. These are available in different 
scales. The details of information available and accuracy 
will depend upon the scale of the map. Such base maps 
provide information of permanent nature e.g. Roads, Tracks, 
Hills and Highgrounds by Contour Lines, Rivers, Canals, 
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Villages and Towns , Forests with type of vegetations etc. 
etc. The military information is then super-imposed upon 
these maps using different symbols. Similarly, the ske- 
tches are made, of the area of interest. However, this 
may show only the military information with bare essen- 
tial ground information. 

1.2. 1.1 Advantages and Disadvantages of Maps 

The advantages of using maps are : 

(i) Maps provide fairly accurate and detailed infor- 
mation of ground. 

(ii) It provides the commander, a 2~dimensional pic- 
torial view in the form of symbols, of the area 
of operation. 

(iii) It is convenient for storage. 

The disadvantages of the maps are : 

(i) The updating is not very convenient. As dis- 
cussed earlier, due to rapidly changing situa- 
tions the deployment has to be changed on these 
maps • 

(ii) The maps are ‘cluttered 1 with unwanted infor- 
mation. The required information has to be 
extracted from it every time it is required. 
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1.2.2 Sand Models 

These are 3 -dimensional extension of maps on a 
limited scale. A sand model is made using sand* cjay* 
wooden and plastic models of various features and mili- 
tary symbols. Sand model provides a 3-dimensional view 
of the battle ground. The tactical military plans can be 
discussed by commanders on these models. 

1.2. 2.1 Advantages and Disadvantages of Sand Models 

The sand models enjoy certain advantages over the 
maps. These are: 

(i) It provides a 3-dimensional view of the ground* 
although on a limited scale. 

(ii) Only such information* which is relevant to the 
discussions* is depicted? thus highlighting the 
features of interest and deleting those which 
are irrelevant. 

The drawbacks of a sand model ares 

(i) The biggest drawback of a sand model is its 

bulk and size. It requires a large quantity of 
stores and a large space on the ground to cons- 
truct it. However this problem is not faced 
where the sand models are not required to be 
shifted* e.g. in training institutions and 
static headquarters. 
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(ii) Its construction is time consuming and labo- 
rious . 

(iii) Because of point at Ser. (ii) above, the modi- 
fication when area of operation changes, is 
also very tedious. 

(iv) Accuracy is limited, and certain features can 
not be depicted to the proper scale. 

1.2*3 Limitation of Maps and Sand Models in Providing 
Non-graphical Information 

So far we have discussed only those features of 
information on maps and sand models which can be depicted 
on them. There is also certain type of information asso- 
ciated with entities on the ground, which can not be 
depicted on the map or sand models, e.g. the load carrying 
capacity of a bridge, population of a village or town, 
unsuitable ground for tank movement etc. etc. We can thus 
partition the entire information system into two parts: 

(i) Graphical information 
(ii) Non-graphical information. 

1.3 Simulation of Sand Model on Computer 

We have seen that present method of using conven- 
tional maps and sand models is not comprehensive and fast 
enough. Whereas the volume of information has increased. 
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the reaction time available with commanders has decreased 
due to fast moving battles of today. Use of computers to 
simulate maps and sand models seems to be the most appro- 
priate solution to this situation. 

Computer-aided cartography has been developed in 
recent years to produce maps. With latest developments 
in hardware and software, soft copy images on CRT have 
replaced printed maps for many applications f ll . Automa- 
tion is also changing the quality as well as amount of 
information that can be provided. A computer aided system 
gives advantages of both, the map and the sand model , 

Some of the features that a computer based sand 
model can provide are : 

(i) Digitising a given map and storing it on a sto- 
rage device of a computer. 

(ii) It can display the map in any selective manner j 

i 

(VL It is possible to display a layer of the | 

i [■' 

map which depicts only the roads or only mili- 
tary symbols. 

(iii) Any desired area can be enlarged to provide 

detailed information. > 

f 

■ 

(iv) Updating, corrections and editing is much simpler | 

i 

and faster. 



e 


(v) The system will provide a central source of all 
the data, which can be selected as per require- 
ment. 

(vi) It can cater for both graphical and non-graphi- 
cal data. 

(vii) Query answering will be faster and more relia- 
ble . 

(viii) Can provide necessary security where only autho- 
rised persons can have access to data. 

(ix) Using animation with raster graphics, movement 
of troops, contingency and different plans can 
be displayed 

(x) Using 'soft copy' display on CRT devices, it is 
not necessary to print the map everytime, which 
saves time. 

1.4 Information Retrieval and Query Answering 

The main purpose of computer-aided sand model is not 
only to create and depict a map, but to retrieve any rele- 
vant information and answer queries. The effectiveness 
of computer-aided maps will depend upon the capability to 
extract information from it, in most efficient manner. 


/ 
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1.4.1 Graphical Data and Queries 

All such data which is useful in the graphical dis- 
play of map on one of the display devices of computer, 
comes in this category. This data will generally be in 
the form of digitised data of entities like roads, rivers, 
contours, and symbols like bridge, temple, survey tree 
etc. This data deal with only 2-dimensional graphical 
queries involving points, straight lines and curves as 
entities. 

Similarly the queries of graphical nature are those 
which make use of such data and perhaps make some computa- 
tions from this data to answer the query. Examples of this 
nature are: 

(i) Distance of a road from a given point. 

(ii) Intervisibility of two points, making canpu** 
tation of heights from contour lines. 

(iii) Width of a river at given crossing point. 

1.4.2 Mon-graphical Data and Queries 

All such data which does not pertain to the display 
in any form comes in this category. This is designed in 
the form of a data base jjjfj. This data base is used to 
answer non-graphical queries, some of which are: 

(i) Whether a given road is under military control 


or civil control 
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(ii) Frequency of trains on a given railway line. 

(iii) Whether the water source is fit for human con- 
sumption or not. 

A query language and data base has been designed 
for this purpose jjfj. 

1,5 Purpose of This Study 

Certain amount of work has already been done towards 
generation of computer-aided map and sand model generation, 
Maj L.K. Chopra , developed a package which can show 
interactively, the deployment of military forces ^3^* 

Capt A.v. Subramanian ' s package generates a base map from 
digitised data £ 2 "]. It also has the facility to create, 
post and delete the symbols, and displays the map at diffe- 
rent layers. 

This study deals with graphical queries pertaining 
to those entities which are represented by curves. It is 
imperative that the system should be fast in making any 
computations, at the same time using minimum of storage 
space. If entire data is stored in the main memory of the 
computer, no doubt it will be very fast, but the amount 
of data, when a set of maps is used, will be very large 
to be stored in the main memory. Therefore proper data 
structure is to be considered. 
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Certain amount of approximation of curves can be made 
whenever the accuracy of query permits «. This greatly redu- 
ces the time complexity of computation* Methods of appro- 
ximation and a hieararchical representation using Strips 
and Strip Tree respectively is implemented. The same has 
been discussed JLn succeeding chapters. Chapter II dis- 
cusses various fceclnniques of approximation of curves by 
breaking them into strips. Chapter III deals with hierar- 
chical representation of strip trees in the form of a 
binary tree, anc3 various operations which can be carried 
out on this structure. Chapter IV gives data structure 
and implem entation details while Chapter V illustrates a 
case study j based "upon the map generated 1 * 1 - 

A list of graphical queries which are discussed in 
this study is given in Appendix A. 



CHAPTER II 


CURVE APPROXIMATION BY STRIPS 


2.1 General 

Two types of information in a map as discussed in 
previous chapter are, graphical information associated 
with the display of map, and non-graphical one, which 
does not pertain to any display but is required to answer 
queries . 

In the remaining chapters we shall discuss those 
graphical entities on the map, which are represented by 
curves. These entities are: 

(i) Roads 

(ii) Cart Tracks 

(iii) Rivers 

(iv) Canals 

(v) Railway Lines 

(vi) Boundaries of specified areas e.g. Forest, 
Marsh Land, Lakes etc. 

(vii) International Boundary and Boundary demarking 
own and enemy forces. 


(viii) Contours. 



13 


2.2 Requirement of Approximation of Curves 

A base map is generated by computer by digitising 
these entities represented by curves. (In rest of the 
chapter a curve will mean one of the entities represented 
by curves.) A digitiser picks up discrete points on a 
curve and stores them as an ordered list of discrete 
points. The number of discrete points picked up will 
depend upon the resolution of digitiser, behaviour of 
curve and accuracy required. Accuracy of representation 
of curve increases with the number of points. 

In order to answer any query pertaining to a curve, 
computations may be necessary. Since curves in a map are 
random in nature and not well behaved regular curves, it 

t 

is difficult to represent them by well defined mathematical 
relations. In order to make any computation, it has to be 
done on the set of points representing it. Considering a 
digitiser with good resolution, and a curve of fairly 
long length, the number of discrete points to represent 
it, may be quite large. Under such circumstances, when 
computations are made on such a curve, considering set of 
all the points, it may consume a large amount of time. 

For example, if we have to find the length c£ a curve 
between two given points, it may involve calculation of 
distance between two successive points and adding them up 
till the other end is reached. Here time complexity of 
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computation is of the order of 0(n) where n is the number 
of discrete digitised points between the two given points 
of the curve.. 

In military applications of maps, it is not always 
necessary to make these computations very accurately, and 
certain amount of error due to approximation may be accep- 
table . This level of error will depend upon the purpose 
for which this computations are made. For example, dis- 
tance between two points, as required for different pur- 
poses can be as follows ; 

(i) In order to site a gun position, and to direct 
artillary shelling on a known enemy position, 
accuracy required may be within a few yards. 

(ii) In order to plan a VHF radio link between these 
two points, an error of few kilometers may be 
acceptable. 

However, this margin of error should be left at the 
discretion of the person using it. 

The strip curves provide a method to approximate 
the curves at different resolutions. 

2.3 Strip Curves 

In order to save computational time, the effort is 
to reduce the number of discrete points on the curve. 
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sacrificing the accuracy to an allowable limit. Concept 
of a strip is to approximate a curve by a straight ’line* 
of thickness w / where w/2 is the maximum allowable error. 
Fig. 2.1 shows such a representation. A strip is a rec- 
tangle enclosing all the points representing the curve. 

The strip with minimum width is called an optimum Strip/ 
which gives minimum error at that level. If the accuracy 
required w r is less than w/2 provided at this level , the 
width of strip can be reduced by spliting the set of 
points (P-^, P 2 ' •••' p n ) into 2 or more disjoints sub- 
sets such that each subset of points results into an 
optimum strip of width w r or less. The methods of such 
sub-division of set of points is discussed in Chapter III. 


2.3.1 Advantages of Strip Curves 

With this approximation/ the points enclosed in 
the strip are virtually reduced to only two. In 
Fig. 2.1/ points (P^ P 2 , .../ P 10 ) are reduced to 
points A and B and/ the curve is equivalent to only two . 
discrete digitised points A and B. Now if the error of 
w/2 is allowed for certain application/ computation from 
these two points A & B will be sufficient/ saving in 
time and storage. (The latter we will see in next 
chapter ) . 
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2*3.2 Notation of a Strip Curve 

In general a strip S can be defined as a 6 -tuple 

01 

(x b- y b' x e- V w l' V 

Fig. 2„'2 shows the general form of a strip. The 
meaning of various terms is as under s 

(Xfo* y b ) - beginning of line B-E, x^rhich is parallel 
to the longer sides of the strip 

(x , y ) - end point of line B-E 

- distance from the B-e to one parallel 
side of strip 

w 2 - distance of line B-E to other parallel 

side of strip. 

(Notation B and E represent points (x^, y^) and 
(x e# y e ) respectively). 

Using this notation/ a curve which is enlosed by , 
this strip can be approximated to the line B-E of thick- 
ness w r « In this case the maximum error for any point 
lying on the curve will be equal to greater of w^ or W£. 

A curve can be represented at resolution w r if there 
exists an ordered sequence of m strips 

k = 0/ ...../ m — 1 
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Such that 

(wj^ + w 2 ) k < w r k = 0, m - 1 

m~l 

X,€ U S 1=1, , n 

k=0 K . 

Where S^. is defined by an appropriate algorithm to 
divide the set of n points / 

X. i=l, ...... n 

into disjoint subsets. These methods are discussed in 
next chapter. 

In above definition, no claim is made of optimi- 
zation of strip to give minumum of error at a given 
resolution. The value of error will depend upon the 
values of w^, and w 2 / which in turn will depend upon the 
method used to form the strip. These methods are dis- 
cussed later in this chapter. 

Another representation of the strip can be made by 
a 5-tuple representation - 

(B, E, w) 

# . 

as shown in Fig. 2.3. In this notation, the line used 
to approximate the curve has been so taken that it is 
equidistant from the two longer parallel sides of the 
rectangle forming- the strip. In this case the maximum 
error will only be w/2. Using this notation, the 
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computations are more convenient, efficient and it ena- 
bles the strip Tree structure with better efficiency, 
as we will see in next chapter. 

2.4 Methods of Constructing a strip 

Dana H. Ballard (jQ has discussed Strip Tree 
structure giving one algorithm for construction of strips. 
Three methods discussed in this study are : 

(i ) Ballard ' s method 

(ii) Least Square distance method 

(iii) Optimum strip method 

2.4.1 Ballard *s Method 

The line approximating the curve is obtained by 
joining first and the last points on the curve. 

(Pig. 2.4). This gives the general direction of the 
curve. The strip is formed by finding smallest rectan- 
gle with two sides parallel to this line, and just 
covering all the points on the curve.. This method uses 
six -tuple notation of strip and maximum error will be 
equal to the distance of the line from the parallel side 
of the strip which is farthest. 

This method, although simple in its implementation, 
gives quite undesirable results in certain odd cases. 

An example of this is shown in Fig. 2.5. Here line B— E, 
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approximating the curve, does not give a very good 
approximation, as also the error will be w^, which is 
very large compared to length of the curve. 

2-4.2 Method of Least Squares 

Erwin Kreyszig has discussed the method of least 
square Distance ^6^ to fit a line through a set of data 
points. 

Consider a set of n points 

<*1' Y l } ' (x 2 ' V' U n' Y n } 

assume that line 
y = a + bx 

is the line passing through the set of points. Vertical 
distance of any point (x^, y^ ) from this line is given 
by (Fig. 2.6) 

d = y . - a - bx j 

Then the sum of square of all the distances is given by 
n 

q = £ - a ~ bxj ) 2 

>1 

The value of q is dependent on values of a and b. For 
value of q to be minimum. 



3a 


and 
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d<Z 

5b 


0 



which gives us the equation of 'Regression Line' as 
y ~ y = b(x - x ) . 

This line passes through point (x, y) with slope b, which 
is called the 'Regression Coefficient* of y on x / *and 
values of x y y and b are given by relations 

1 

X SB — (x, + X ) 

n 1 
1 

y - ~ ( n + 

n 


Sxy 


S 


Sxy is called 1 Convariance of Sample * . Here 


n 


n~l 


Z<*j - 


x) 2 


J“1 


n 


and Sxy = 


— 5 1 (x. 

n—l ^ 

n 1 j=l 


xHyj 


- x) 


Therefore by computing values of 3^, Sxy we can 
compute slope of the line which passes through (x, y). 
The strip can then be formed as in Ballard's Method by 




finding two sides parallel to this line and just enclo- 
sing all the points. This also uses' the 6-tuple notation 
of the strip. It was experimently seen that strips such 
formed are not optimum strips. 

2.4.3 Optimum Strip Method 

An optimum strip has been defined as that strip 
which has a minimum width for given set of points defi- 
ning the curve. By ensuring that the width is minimum, 
the margin of error is reduced to minimum, The method 
uses the notation of 5-tuple representation for the 
curve. Fig. 2.7 shows an optimum strip constructed by 
this method. 

Like previous method of Least Square Distance, the 
dsape of curve has no meaning, only the set of points 
describing the curve are considered. * 

This method consists of two parts: 

(i) To determine the direction of optimum rectangle 
(ii) To construct the strip. 

2. 4. 3.1 To Fjnd Direction of Optimum Strip 

This involves computing smallest width for the 
given set of points, at different angles and comparing 
the results. 

Consider a set of points X^, ...... X n descri- 

bing a curve, and a line 





ax + by + c = 0 


If we compute perpendicular distance of each point from 
this line, then the two points which are maximum and 
minimum distance from this line respectively, mark the 
two boundaries of the strip. The difference of these 
distances d^^ and give the width of strip which 

is parallel to this line, the two sides passing through 

P dmax and P dmin res P ect: *- ve ly • Fig* 2.8 illustrates this 
method. In Fig. 2.9, point Pg is at max distance from 
given line, and point at min distance. (Only magni- 
tude is not sufficient for this, thus Pg, and P^q which 
have lesser magnitude of distance but with +ve sign are 
at greater distance than that of P^). Therefore, the 
two dotted lines mark the two parallel sides of the strip 
which is constructed in the direction of given line. 

The next step in this part involves repet ing this itera- 
tion for lines at different inclinations. The algorithms 
for choslng steps and different values of directions is 
discussed in next chapter. 

For the ease of computations, we consider the 
lines which pass through the origin of the coordinate 
system. Thus we have lines 

y = mx 

for different values of m. Distance of a point (x^, y-^) 
from the line 



ax + tox + c = 0 


is given by 

axi + by^ + o 

d= 

but if the line is y = mx, passing through the origin, 
then 


a = m ■» Tan 9 
b = -1 

(refer to Fig. 2.7). Then the distance is given by 
d « x 1 sin 9 - y^. cos 9 

By changing the value of 9 from 0° to 180°, in the 
steps as determined by suitable algorithm, we can compare 
the width of strip in each direction. The strip having 
minimum width will give us the direction of optimum 
st rip , um • 


2. 4. 3. 2 Constructing Optimum Strip in 


Direction 


Having found the direction of Optimum strip, we 
get two points through which pass the two sides of rectan 
gle which are parallel to the direction. The other two 
sides of the rectangle can similarly be found by finding 
the two points which are at maximum and minimum distance 
respectively from a line which is perpendicular to the 


direction of strip. Line B— E approximating the curve 

is obtained by joining the midpoints of these two lines 

which are obtained in the direction perpendicular to the 

optimum strip direction. Fig. 2.7 shows an optimum 

strip with four boundary points , and P^ lying on the 

sides which are parallel to the direction. P and P„ on 

2 4 

the sides perpendicular. Line B-E will approximate any 
curve which has points lying within this rectangle and 
includes P 1# P 2 » P 3 and P . Therefore we observe that 
ordering of points, or in other words shape of curve 
has no effect on the strip. 

2*4.4 Comparison of the Three Methods 

The basic purpose of using strips is to approxi- 
mate a curve such that it reduces the computation and 
storage effort, at the same time with minimum possible 
error. In Ballard Method, the error in worst cases, may 
be very large as illustrated in Fig. 2.5. This does not 

t 

give a true approximation of the curve. 

The Least Square Distance method is improvement in 
Ballard Method where shape of curve is neglected but only 
the set of points on the curve are considered. However, 
it involves larger amount of computations, and may still 
not give an optimum strip. 

The optimum Strip Method also involves larger 
amount of computations, computing distance, of all the 
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points from lines at different inclination,, but it gives 
the best approximation to the curve, -although there 
might still be a quantisation error due to steps of angle 
from 0° to 180° chosen for various lines , but this error 
will be negligible if these steps are chosen properly. 

In this study t he method of Optimum Strip has been 
used to construct the strips. It's implementation details 
are discussed in succeeding chapters. 
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CHAPTER III 

STRIP TREES AND BASIC OPERATIONS ON STRIP TREES 

3.1 General 

In Chapter II, methods of approximating a curve by 
a strip were discussed. The strip enclosing the entire 
curve is called the 'Root Strip*. The width of the Root 
Strip will depend upon the set of points defining the 
curve. In case of Ballard's method of forming a strip 
it also depends upon the position of first and the last 
point in the ordered set of points. This Root Strip is 
not of much practical use since the width of this strip 
is not predicatable , and hence the level of error is 
also not predicatable. In order to reduce the width of 
the strip to an acceptable error level, we have to split 
the set of points to smaller disjoint subsets. 

If 

S = ^p 1# P 2 , P n ] 

be the set of points defining a curve, then we can parti- 
tion these points into smaller subsets 
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Fig 3.1: Representation of Curve by Strips 
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such that 


Sl “i P l' P 2 P i 1 

®2 = [ p i- *>1+1 Pj] 

} 


p n] 

and 


w sk^ w r 


where 



w sk is width of strip for Subset 

w r minimum required width for allowable error in 
approximat ion . 

Figure 3.1 shows an example where a curve defined , 
by 15 points has been subdivided into 4 strips, each 
strip having a width less than or equal to w r the requi- 
red minimum width. The curve, which was represented by 
15 points can now be represented by these four strips. 

If still lower resolution, say w^/2 is required, then 
first, second and third strips will be again subdivided 
into strips of lesser width. 
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3.2 Strip Trees 

A strip Tree is a hierarchical representation of 
a curve in the form of strips, such that each lower 
level corresponds to finer resolution representation of 
the curve. This is achieved by considering each subset 
of points as a new set, and subdividing each set of 
points whose Hoot Strip (also called Hoot Node) has 
width greater than w . This tree structure has proper- 
ties desirable for various computations required to be 
carried out on these curves. Intersection and point 
membership (for closed curves) can be resolved with time 
complexity O(log n), where n is the number of points 
describing the curve. The curves can be efficiently 
encoded and displayed at various levels of resolution. 
The representation is closed under intersection and 
union and these too, can be carried out at different 
resolutions. All these properties depend on the hierar- 
chical tree structure which allows primitive operations- 
to be performed at lowest possible resolution with great 
computational time savings. 

3.3 Methods of Forming a Strip Tree 

Efficiency of Strip Tree structure for computation 
will depend upon the method of representation of the 
tree. Theoretically, each strip with width w r can be 
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subdivided in n— subtrees, thus forming a n-ary tree struc- 
ture. Basic consideration in methods of sub-division is 
that leaf strips of width w r should be formed as fast as 
possible, keeping in view the computational complexity 
in mind. The four methods of sub-dividing a strip are: 

(i) Tertiary Strip Tree 

( ii ) Maxpoint Binary Strip Tree 

(iii) Midpoint Binary Strip Tree 

(iv) Leaf -Node- First Strip Tree 

3.3.1 Tertiary Tree Method 

Fig. 3.2 illustrates this method. We have seen in 
Chapter IX, that while forming a strip, each of the four 
sides of rectangle will pass through atleast one point 
each on the curve. This is to minimise the size of strip, 
after the direction of B-E line has been found. Tertiary 
Strip Tree Method uses two points which lie on the lines 
parallel to B-E line, to sub-divide the strip. 

If strip S is enclosing a curve with points 

S = ^P ± , P 2 , P n} 

then its sibiaigs will be the strips enclosing points 

S 1 * i P l P maxl^ 

s 2“i P maxl * P max2^ 

s 3 “ ^ P max2 P n] 










where p max i and p maX 2 are t * ie two P° : >- nts on the li nes 
parallel to B-E and p maxl -is earlier than P max2 in orde- 
ring of points P^.... P . 

Tne method makes use of probability that the curve 
has a more uniform slope in these three parts than in any 
other position. Convergence to the leaf strips of width 
w is fairly fast in this method. However , a tertiary 
tree is more combursome for making computations, although 
more efficient to form. 

3.3.2 Maxpoint Binary Strip Tree 

This method is a Binary-Tree adoption of previous 
method. In order to reduce the tertiary tree to a 
binary tree, for ease of computations, it uses one of 
the two points on the lines parallel to B-E (Fig. 3.3). 
Here a consideration, as to which of the two points is 
to be selected as dividing points, has to be made. 

In case p maxl is such that ifc comes tefore p m ax2 '' 
in the order of points, and it is selected, then strip 

for points 

s 1 * 1 P 1 P maxl ^ * 

will converge faster because it will have lesser width 

than other one. 
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<P m«l , I & (p max2 ) I “ Siting points in first 

sub-division 

^ P maxl^II & ^ P max2^il “ Limit:Ln 9 points in second 

sub-division 


\ 

\ 

r 

I 

i 


| 

| 


Pig 3.4: Illustrates Inefficiency of Max-point 

Strip Tree Method , where P max ^ is chosen 

as Dividing Point without Considering 
its Location. 
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Fig 3.5: Mid-point Sub-division Strip Tree j 

Representation s 
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This method can become very inefficient in the 

worst case as illustrated in Fig. 3 . 4 . Here in first 

two iterations, the left subtree has only 2 points, 

because in each subtree second point happens to be P , . 

“maxi 

To avoid this situation, a test can be made on P 

maxi 

and p max 2 * The points which lies nearer to the line, 
perpendicular to B-e and passing through the mid point 
of B~e should be selected as subdivision point* This 
will ensure that the point is not close to either first 
or last point of the curve lying inside the strip. This 
method has to make extra computations and comparisons 
thus making it slower. 

3,3,3 M *<3~PCdnt Binary strip Tree 

The method is simplest and therefore fastest in 
implementation. If the strip has n points then, sub- 
division point p^ is given by 

»d* p ,ru neven 

Pa'V, n 043 ■ • 


Therefore both the sub-strips have either equal number 
of points or a difference of one point. Complexity of 
this method will be Q(log n). (Fig. 3.5) 



(c) 


Fig 3.6s 


Leaf-First Strip 
(a) Leaf Nodes 
(c) Root Node 


Tree Method 
(b) Intermediate Step 


* 
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3,3*4 Leaf -First Strip Tree 

50 far we have discussed only top down approach 
of forming a strip free, where, each strip yields two 
siblings. This method is reverse of top-down technique, 
in which two adjoining strips are combined to creat 
their parent. 

Initially, the algorithm starts with constructing 
leaf nodes, each having width 4w r# This can be done by 
selecting subsets 

5 1 - l p l* P 2 p 4 

5 2 “ i P i P j ] 

s 3 * i p j ^ 

s p ■ i p n 3 

such that 

w(s k ) i w r 

and w<s )c , first point of s k+ i^> w r < k = 1, ..... p, where 
w(« k > is width of atrip. £orj*ad for points Of 

Having formed the leaft nodes, we can form nodes 
with coarser resolution. by combining points of two succe. 
ssive nodes and forming a single new node. This itera- 
t±on stops when a single root node is formed. Fig. 3.6 
illustrates this, where initially 4 nodes were formed. 



aha in twp, successive iterations 2 and single root node 
is formed » 

3.4 Basic operations on strip Trees 

Computations on a curve can new be performed using 
the strip Tree formed for the same* The computational 
complexity of various operations which can be performed/ 
may be difficult to characterize because it depends on 
the particular geometry of the curve. If the curve is 
well behaved i.e. if it is relatively smooth/ then the 
algorithms are very efficient. If we are using binary 
trees/ and look at only one of the two sons at any node, 
then the time complexity is of the order of O(log n). 

The curves generally encountered in a map fall in this 
category. There are some closed curves, however, which 
are not well behaved e.g. contours, boundary of speci- 
fied area etc. But in such cases, the nature of queries 
generally are of point containment which, as we will 
see later, can be computed quite efficiently by using 
higher level nodes itself. 

The operations which can be performed on Strip 
Trees, and are used in query answering techniques of 
Chapter IV are given beiow. This assumes that a binary 
tree with leaf nodes of width w r has been created. 
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3.4.1 Display of Curve at Different Resolutions 

As discussed earlier, a curve may be represented as 

a set of strips such that each strip has a resolution 

width less than some fixed value w . If we have a dis- 

r 

play device depended procedure STRIPDISPLAY to draw, a 
rectangle, then algorithm to display a curve at resolution 
will be as given below. (A pseudo~Pascal language is 
used to describe all the algorithms.) 

procedure DISPLAY (T, w ); 
begin 

if w(T ) < w then STRIPDISPLAY(T) 
else DIS PLAY (LSON(T ) , w r ) and 
DI SPLAY ( RSON ( T ) , w^J 

end; 

3.4.2 The Length of the Curve 

The Strip Tree provides a very simple method of 
calculating the length of a curve at given resolution w^. 

function LENGTH (w r , T) sreal; 
begin 

if w(T)^w r than 

LENGTH :® (SQRT((s b (T) -x e (T)) + (y b (T) 

y e CT)) 2 ). 


else 
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Pig 3*7: Algorithn for Computing Length using 

Strip Length Causes Unpredictable Error 




! 
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LENGTH := LENGTH (w , LSON(T)) + LENGTH (w , 

* X* 

RSON(T) ) ) 


end; 


This algorithm sums up the lengths of all those 
strips which are just equal to or less than w , while 
travelling down the tree. This may cause unpredicable 
error in case of such curves which have too many short 
turns, as shown in Fig, 3.7. Once again the assertion 
is made, that if the curve is well behaved, the approxi- 
mation may be acceptable. However in actual implementa- 
tion, a modified algorithm where actual length of curve 
is stored in each strip, has been used. The same will 
be discussed in next chapter. 


3.4.3 Intersection of Two Curves 

One of the important features of the representation 
is its ability to compute intersection between curves 
at higher resolution itself. In order to explain the 
algorithms for intersection, following definitions are 
necessary. 

(a) Two strip segments, S 1 derived from curve 

C. and S n derived from C 0 have a ‘null' inter- 

section iff s iH S 2 ~ 

(b) Two strip segments and have a ' clear * 
intersection iff all the sides of the strip 
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parallel to line segment (B-*E)^ and (B-E)^ 
intersect . 

(c) Two strip segments and S ^ have a ' possible 1 
intersection if condition (b) is not satisfied 
and yet S 2 # 0. 

These cases have been illustrated in Pig. 3.8. A 
fairly obvious but very important lemma is givens 

CLEAR INTERSECTION LEMMA: If two strip segments have a 
clear intersection and the strips are both regular# then 
the corresponding curves must also intersect. Fig. 3.9 
shows such an intersection. 

From this lemma we can make the assertion that if 
two root strip segments corresponding to two curves# do 
not intersect then the curves do not intersect# and if 
the root strip segments have a clear intersection then 
the curve also intersects. This important assertion 
helps in finding intersection at higher levels itself# , 
without having the need to travel the trees to leaf 
nodes of width w r or less. 

in the algorithm given below, function STRIP- 
INTERSECTION determines the type of intersection between 
two strip segments. A(T) is the area of strip at node 

T. 



m 


Closed Area 


Fig 3.10: To test whether a Point is inside a 
Closed Area 


Closed Area 


/ Closed 
/ Area 


3.11: Cases to illustrate i-urihaguataes an Testang 
Containments a Point in a Closed Area . 
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function INTERSECTION (T^ T 2 ) s boolean; 
begin 

Type SB STRlPINTERS£CTION(T., , T ); 

Case Type of 

Null; INTERSECTION := false; 

Clear: INTERSECTION := true; 

Possible: if A(T 1 )>a(T 2 ) then 

INTERSECTION := ( INTERSECTION (LSON(T' 1 ) / T 2 ) 

OR ( INTERSECTION (RSON(T 2 ), T 2 ) 

else 

INTERSECTION := (INTERSECTION^ , LSON(T 2 )) 
OR (INTERSECTION^, RSON(T 2 ) ) 
end; ^case ^ 
end; 


3.4.4 Determining whether a Point is inside a Closed 


Curve 


The Strip Tree representation of an area by its* 
boundary allows in a simple and straight forward manner, 
determination whether a points is inside a closed curve. 
Following Lemma is used in the algorithm for determining 
the containment of a point 1 



POINT CONTAINMENT LEMMA: If any saini-infinite line termi- 
nating at the point intersects the boundary of an. area 
an odd number of times, the point is inside the area. 
(Fig. 3.10) 
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To decide whether a point P is member of an area 
represented by a Strip Tree, we need only to compute 
the number of clear intersection of the Strip Tree with 
any semi-infinite strip, S Q which has width 0 and emana- 
tes from P. If this number is odd then the point is 
inside the area. The algorithm given below uses a proce- 
dure CREATSTRIP (S^, P) which creats a strip for semi- 
infinite line emanating from P. 

function INSIDE (p, T): boolean; 
begin 

CREATSTRIP ( S Q , P); 

NUMBER := NOINTERSECTIONS (S Q , T); 

if NUMBER is odd then INSIDE := True 
else INSIDE := flase 

end; 

function NOINTERSECTION (S, T) : integer; - 
begin 

TYPE := STRIPINTERSECTION (S, T); 

Case Type of 

NULL: NOINTERSECTION := ); 

CLEAR: NOINTERSECTION : = IM- 
POSSIBLE: NOINTERSECTION := NOINTERSECTION 
(S, LSON(T)) + NOINTERSECT ION(S, RSON(T)) 

end; % case \ 


end/ 


54 


This algorithm runs into trouble when strip S Q 
is tangent to the curve (Fig. 3.11(a)). Similarly if 
strip S Q passes through an end point of a node in strip 
Tree corresponding the curve , number of clear intersec- 
tion will be 0/ causing an ambiguity. (Fig. 3.11(b)). 
These ambiguities can be, however, overcome by more 
complex algorithms. 

3.4.5 Other Operations 

Besides the operations discussed so far, certain 
other operations which can be performed on Strip Trees 
are : 

(i) Intersection of a curve with an Area 
(Ii) Intersection of Two Areas 
(iii) The union operations. 

(i) and (ii) are essentially versions of intersec- 
tion of two curves where Areas are described by closed 
curves. The operation of union has not been used for any 
of the query answering algorithms of Chapter IV, and 
therefore it has not been discussed here. 
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CHAPTER IV 

ALGORITHMS, DATA STRUCTURE AND IMPLEMENTATION 

4.1 General 

In Chapter III, techniques of forming Strip Trees 
and certain basic operations which can be performed on 
Strip Trees, were discussed. In this chapter we shall 
discuss the design of a package and its implementation 
for processing various graphical queries making use of 
Strip Trees. The scope of queries is restricted to 
Military applications only. 

A compromise between efficiency in time and sto- 
rage space is essential for any such system, and a suita- 
ble Data Structure is necessary for the same. Dynamic 
storage allocation, using linked list, has been used 
wherever possible. To make the package, user friendly, 
a suitable command language has been designed. 

4.2 Scone and Capabilities of. Package 

The package is developed to process all such 
graphical queries of Military nature, which can be 
answered using existing graphical data on a computer-aided 
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map. This data is essentially in the form of curves 
describing various features on the map, as discussed in 
earlier chapters. 

Tho queries which can be processed are: 

(i) Distance. 

(ii) Length of feature. 

(iii) Containment of a point within a closed area. 

(iv) Area of closed areas. 

(v) Height of a point on the map. 

(vi) Intervisibility between two points. 

(vii) Magnetic, Grid and True bearings from a point 
to another point. 

(viii) Intersection of features. 

The various modes of these queries are possible. 

A list of such modes and corresponding user commands, is 
listed at Appendix 1 B '. 

The package has been designed to be compatible 
with map generation package earlier developed \ 2 ~\. 
uses the same input source for creating strip Trees for 
various features, which has been used to generate the 
map. This input is in the form of digitised data for 
curves, available in source file DATaIN. A serial number 
has been allocated to each feature within the same class 
of features. This enables package to distinguish from 
one feature to another. 
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The command language gives the flexibility to user 
to give arguments to query in any order (except in cer- 
tain cases). An error detection system, embedded in 
the package, can report the error to user at different 
stages. Package accepts all commands in strings form, 
to avoid error in data type in input. The error in data 

type, if any, is detected by the program and appropriate 
error message given. 

4 * 2 * 1 Ref erence Technique of Referencing a Point 

A six figure Grid Reference (gr) technique of 
referencing a point on the map, which is commonly used 
in Military usage, has been used for interaction between 
the user and the system. GR is a 6-digit integer giving 
location of a point on the map. A GR 457695 represents 
a point on the map with coordinates (45.7, 69.5) using 
the Grid coordinate system as depicted on the map. 
Whenever a user has to specify a point or the system 
gives a reply to a query in which a point is involved, 
6-figure GR is used. 

4*3 Organisation of Package 

The package is divided in 3 functional modules. 
These modules are independent of each other and perform 
specific tasks. Module CREATETREE constructs strip Tree 
for a feature (or set of features) for specified w r * 
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Module COMMANDPROCESSOR consists of Command Language 
Processor (CLP), it analyses the command and passes 
control to appropriate part of QUERYPROCESSOR module. 
Finally, QUERYPROCESSOR module performs all the necessary 
computations according to the nature of query and gives 
results. The modular layout of system provided flexi- 
bility to modify any module without affecting others, 
betails of implementation of each is described in succee- 
ding paragraphs. Fig, 4.1 shows the functional organi- 
sation of the package. 

4.4 Module CREaTETREE 

This module takes digitised data of a feature as 
the input and creates corresponding strip Tree, in the 
form of a binary tree. The leaf node strips of binary 
tree have width ^ w r , which can be specified by the user. 
This module is in the form of an independent program and 
makes use of the fact that creation of Strip Trees is a 
one time process. In other words. Strip Trees need to 
be constructed only once before making use of other two 
modules. Thereafter, this module is not required until 
a Strip Tree corresponding to some additional feature is 
to appended, or the resolution of existing ones is 
required to be modified. 
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Basis of chosing algorithms and data structure is 
guided by the fact that it is only a 'one time use* 
module. Accuracy, rather than efficiency in time, is 
of more importance, since all further query processing 
is to be carried out on Tree structure. 

The module uses 1 optimum Strip Method* for making 
strips at each node and 'Mid Point Binary Tree* for 
Strip Tree representation. 

4.4,1 Data Structure 

The computations in this module are performed on 
the set of points describing the feature, strips formed 
at each node have all the necessary information which 
is required for subsequent computations. 

A buffer PT, which is a 2-dimensional array of 
[1000x2] , i« used to store x and y coordinates of 
points describing the curve* These points are obtained 
fran file DATAIN, for one feature at a time, and stored > 
in this buffer. PT[l,x] and pt[i^y] give x and y coor- 

4" Vi 

dinates respectively of the 1 point on the feature. 
This indexing of points is useful while dividing the 
set (or subset) of points in two halves as required for 
Mid Point Binary Tree Method. 

Internal representation of strip Tree is done by 
using a linked list* Each node of the tree contains 



4. 4.1*1 Freelist 


v t. 
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The program is capable of processing data for any 
number of features. In such a case, unless the nodes 
used for each feature are not released for the subsequent 
use, they will keep getting accumulated. The size of 
storage space required will be 

n 

0 ( Z up 

i=l 

where is the number of nodes created for i * feature 
at given resolution, and n is number of features. This 
size will increase with increase in .number of features. 

Maintaining a FREELIST of nodes, overcomes this 
problem. After the Strip Tree for a feature has been 
built, it is stored onto file STREE and all the nodes 
used in the TREE are released to join a list, called 
FREELIST. Whenever a node is required for adding a node 
to Strip Tree (except in case of ROOT node of a tree, 
it will always be a pair of nodes corresponding to left 
and right sons), first a check is made in FREELIST. In 
case the list is not empty, node is taken from FREELIST, 
otherwise it is asked from system, using new (P). The 
order of space required for nodes is 0 (max ^ ) . This 
is not dependent upon number of features. 
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4*4,2 Building a Strip Tree 

Optimum Strip Method as discussed in Chapter II 
is used for building a strip for a given set of points. 
The set of points is specified by the indices of first 
and last point in array pt. The details of some of the 
important procedures/functions in this module are given 
below. 

(a ) Procedure DIRECTION ( INC :real ; FIRST , LAST ;integer ) ; 

This procedure finds the direction of optimum 
strip for points which have indices from FIRST to LAST 
in Buffer array / Pr. The procedure computes width of 
strip for these points at different angles between 
MINANGLE and MAXANGLE (these are global variables ), at 
an incremental step given by INC. This procedure is 
called twice for every set of points. First time it is 
called with values of MINANGLE and MAXANGLE as 0.0° and 
180.0° respectively, and value of INC as 10.0°. If 0 
is the direction in which strip has minimum width in 
first iteration then in second iteration, value of 
MINANGLE and MAXANGLE is set as 6 + 15.0° and INC as 1.0°, 
and process repeated. This gives an optimum strip within 
an error of + 1.0°. 
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(b) — Ction FWlm,H real /FIRST, LAST; Integer): real ; 

This function is called by procedure DIRECTION to 
compute the width of strip for a given set of points 
(indices FIRST to last) and an angle THETA, it computes 
distance of all the points from line - 

y « mx 

whore m = tan (THETA). 

<c) Procedure MAKESTRIP ( FIRST , LAST : integer ) ; 

This procedure calls procedure DIRECTION with two 
sets of values of M INANGLE, MAXANGLE and INC as discussed 
earlier and builds an optimum strip in direction computed 
in second iteration of DIRECTION. It also computes four 
corner points of the strip. 

(d) Procedure CREATRBB ( FATH m : PTR; FIRST .LAST ; integer ) ; 

It is a recursive procedure to build the Strip 
Tree. Each time, it makes a strip for points with indi~ 
ces FIRST to LAST by making call on MAKESTRIP. If the 
width of this strip is w then it makes a recursive call 
on itself to create left and rightsons as - 

CREATREE (LSON, FIRST, I) 

and CREATREE (RSON, I, LAST) 

I » (FIRST + LAST) DIV 2. 


where 
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In this manner it creates a depth-first bin ar y 
Strip Tree . 

(«) Procedure primtmode (p?ptr); 

This is also a recursive procedure which outputs 
the Strip Tree nodes onto file STREE. The data in file 
STKEE is stored with first line as feature identification/ 
followed by nodes. Identification line has the format: 








■ ■ 

Feature 
Ser. No. 

NODES 

COUNT 

MINWR 

MAXWR 


Where 


Feature Class - Class or feature e.g. RIVER. 


Feature ser. No. - Ser. number of feature in its 

class e.g. RIVER 5. 

NODES - Total number of nodes / including 

Root node and leaf nodes, in the 
Tree. 


COUNT 


MINWR 

MAXWR 


- Number of digitised points of 
the curve. 

- w . for which the tree is built. 

r 

- Width of Root node. 
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Subsequent lines (number of such lines is equal to 
NODES) contain information about each node. Their 
format is : 


dER * NO* 

OTDTH 

LENGTH 

LS ON SER 

CORNERS 


(f ) ether i-rocedures 

Certain other procedures are as under : 

- hibONaMb: identifies the name of feature from 

data file DATA IN. 

- INjaTAs reads digitised data from DATalN and 

stores in buffer array PT. 

- GDTNdJ: Gets a node either from FREELIST, if 

available, or from system. 

~ FREELIST: maintains a list of free nodes. 


4.5 Module COMMANPPKOCESSOR 

As discussed earlier, the package has a user 
friendly command language (CL). The purpose of the 
language is to make the query commands as simple as 
possible, and at the same time powerful enough, to avoid 

any ambiguities. 

The package is designed, keeping in view the fact, 
that users who are going to make use of it, may not 



67 


have ]cnowledge of any programming language. The user 
should not be burdened with any language specifications 
or system depended features. The entire package , the 
atrip Tree Representation of Curves , the degch of 
resolution etc. should be transparent to the user. 

4.5.1 Features of CL and CLP 

(i) Comm ant., s are simple , with a syntax quite simi- 
lar to English. For instance, for a query: 

'find dISTaNCE from point 234569 to a feature 
KO.'PriH NO. 5#' 
corresponding command is 

.. ICT-'U'ICiV , 2 34 56 9 / RQaNNH. 5 

(ii) abbreviations of feature class and query name 

are permitted, a feature class can be abbre- t 

of* f 

viated to its first 5 characters, whereas in \ 

case of query name it can be abbreviated to 

its first ^ character. For example, CaRTTR jj 

, i 

for CdlTTRACK, CIST for DISTANCE can be used. I 

I 

\ 

(iii) Referencing of points is done using 6 digit | 

* 

Grid Reference (GR), commonly usee,, in Mill- j 

fe 

I 

tary applications. i 

(iv) It provides flexibility in the ordering of j 

i 

arguments to query . For example •: 
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DIST/2468 92/KMER . 5 
and 

DIST/RIVER.5/246892 
will be interpreted as same query. 

(There are certain exceptions to this# which are 
discussed later). 

(v) 1 HELP* Command 

This command provides necessary help to the user 
about usage of this package, whereas HELP pro- 
vides help of general nature# HELP.eactn 
where extn is the query name# provides help, on 
specific type of query specified by the extension . 

(vi) *7* Mode 

Once the exeuction of program begins# the program 
always remains in this mode. A' 7 prompt on TTY 
indicates that program awaits user command* User 
command can either- be a HELP command or a QUERY 
command. After every command# the program proce- 
sses it# and after one of the following two condi- 
tions# returns to this mode - 

(a) Giving HELP message dr reply to the query# as 
the case fttey * ojt 

(b) Giving an error message# if there is an error 
at any stage. 
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(vii) 1 Exit 1 Command 

This command in *?* mode terminates the exe- 
cution of . pqrogram . 

(viii) Blank Suppression 

CLP suppresses all blanks in user commands, . 

(ix) Upper Case - Lower Case Alphabets 

User can use upper or lower case alphabets. 
However, CLP converts all lower case letters 
to upper case and all processing is carried 
out on upper case alphabets, 

(x; ) User commands are accepted as character strings 
by .CLP, and converted to integer or real num- 
ber wherever necessary. This eliminates any 
errors caused by conflicts of data -type in 
the input, which if any is detected by CLP and 
reported* 

4 .5.2 Crammer of CL 

BNF grammer of command language (CL) is given at 
Appendix ' D 1 » 

4.5.3 implementation 

As discussed earlier , CL accepts the user command 
input as a character string; suppressing all the blanks. 
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A BUFFER of 80 characters is used to store the user com - * 
mand read from TTY, after supressing all the blanks. 

Further processing is done on data taken from this BUFFER. 

*/' is used as delimeter between query name and arguments. 

The nature of query is then obtained from the BUFFER* 

Following important procedures/functions are used upto 
this stage ; 

- procedure READLINE : Reads command from TTY and 

stores in BUFFER. 

- procedure TESTQUERY :Tests for nature of query/ 

HELP command. In case of 
■ HELP or help. E xt n, gives app- 
ropriate Help message. In 
case of a query, gets a serial 
number corresponding to the 
query name* 

f 

- function FINDQN: Returns query number as an 

integer corresponding to the ' j 

serial number of the query in j 

| 

the ordered list of query names, [ 

f 

as given at Appendix *E*. j 

\ 

4 *5.3.1 Code Generation 

Cnee CLP has obtained the information regarding 
the nature of query (hot a HELP command ) , it has to 
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obtain its arguments and process them for their validity 
etc* The arguments which follow the query name can be 
of following types: 

(a) GR.of A Saint ; 

It is a 6 -digit integer number corresponding to loca 
tion of a point. 

(b) Feature Name 

It has name of feature-class , followed by the Serial 
number of feature under that class, e.g. in case of 
RIVER No. 7, RIVER is the name of feature class to 
which It belongs and, 5 its serial number under the 
class. This is specified as RIVER, 5 in the command. 
When any or all of the features under a feature- 
class are to be referred, ANY or ALL is given in 
place of serial nunber (e.g. RIVER. ANY or CARTTRACK. 
ALL) . The feature can further be subdivided Into 
2 groups* 

(i) Open Features - River, Road etc. 

(li) Close Features - Contours, Boundaries of areas 

etc . 

This distinction is made after an op-code is gene- 
rated cor re sponding to the query (explained later). 
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(c) Special /'arguments 

Certain special arguments which pertain to specific 
query or argument can be given. They are: 

- INSIDE: To specify portion of a feature 

lying inside a closed feature. 

- ENTIRE: In LENGTH query# when the length of 

entire feature is required. 

« MAG/ GRID/TRUE: To specify the type of Bearing to 

be computed. 

A 7-digit integer op-code is then generated for 
uaC h argument. The code has following information: 


1 2 4b L- 

TYPE 

CLASS NAME 

SER. NO. 


Digit 1 - first digit of op-code gives type of 
argument# as given below : 

0 : a 6-digit integer GR 

1 : an open feature 

2 : a closed feature 

3 • special argument 

Digit 2 to 4 - Digits 2# 3 and 4 give serial num- 
ber of feature-class in the ordered list 
of feature-class# or the serial number 
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of special argument from the ordered list 
of special arguments , as given at Appen- 
dix 'E*. In case of Type GR, these digits 
combined with digit 5-7 give 6 -digit GR. 

Digits 5 to 7 - In case of a Type Open or closed 

feature , digits 5-7 contain the Ser. No. 
of feature in its class. In case of ANY 
or ALL# value of these digits in 000 and 
999 respectively. For Type Special argu- 
ments/ value is 000, and in case of Type G 
GR, these digits, combined with digits 
2-4 give 6-digit GR of point. 


Examples i 



Aroument 

Op-Code 

(i) 

Point 245659 

- 0245659 

(ii) 

river. s 

- 1004005 

(iii) 

RIVER. ANY 

- 1004000 

Uv) 

C0UNT0UR .ALL 

- 2001999 

(v) 

INSIDE 

. - 3002000 

(vi) 

HAG 

- 3003000 


4 . 5. 3. 2 Analysis of Query, 

This part of. module COMMAND-PROCESSOU, carries out 
on analysis of argents along with the nature of guery 
asked for. Xt performs following takes 
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(1) Tests for correct number of arguments as per 
the type of query. 

(i±) Tests whether arguments are compatible with 

type of query. For instance , in case of query 
'HEIGHT', argument should only be a point GR. 


(iii) In case of argument as feature name# it tests 
whether feature should be open or closed fea- 
ture . 

(Iv) Tests whether ordering of arguments# where 

flexibility of ordering is not permitted, is 

correct . 

nfter testing various conditions, it re-orders the 
arguments where the flexibility to the order is allowed, 
and the same has been changed. It then makes call on one 
of the 'Buffer* procedures with appropriate arguments and 
passes control to module QUERY-PROCESSOR. 


4.6 Module Q UERY-PROCESSOR 

This is the part of package which performs all the 
necessary computations to answer the query. It is initia- 
ted by one of the -Buffer' procedures called by COH.-MD- 
PROCESSOR module. Each of these buffer procedure corres- 
ponds to a specific type of query as specified by set of 
commands in Appendix 'B'. The module makes 
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In Strip Tree form in file STREE, and creates an inter- 
nal representation of Strip Tree or list of leaf nodes 
upto those nodes of strip Tree in STREE which are just 
specified by user for computations. 


4-6.1 Data Structure 


As in the case module CREATTREE, a linked list 
using pointers, has been used for internal representa- 
tion of the TREE or LIST of the nodes required for the 
computations!, Oily one array - 


CHECKTABLE j packed array 


[1...500] 


of boolean? 


is used, whose purpose is explained later. 


4-6,2 implementation 

As discussed earlier, depending upon the nature 
of query, either a LIST of those nodes in the tree which 
are just < w , or an internal representation of the Strip 
Tree upto the resolutipn w r is required to be built. 

The purpose and algorithm to build a LIST or TREE is 
explained# 

4,6, 2.1 LIST of Nodes 

For the queries, DIST, LENG, CROS (where point of 
intersection 1 a also required), AREA, HEIG, INTE and CON- 
TAIN, a LIST of those nodes in the Strip Tree, which 
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contains strips of width <w , is required. This is 

Ni' 4. 

maintained as a Linked List. 

Algorithm used, performs a sequential search of 
nodes starting from the ROOT node. Whenever a strip 
which has the width w is encountered, it is- made to 
join the LIST. CHECKTABLE is used to ensure that, those 
strips which are SONS of a node joining the LIST, and have 
width ^w^., are not included, when encountered subsequently. 
This is achieved by making entry corresponding to such 
SONS as 'false* in the CHECKTABLE. 

4 .6.2.2 Creation of TREE upto Resolution 'w 1 

Certain queries require the tree structure of the 
strips for efficient computations. These are CONTAINMENT 
and CROSSING. In order to have, an internal representation 
of Strip Tree using linked list, an auxiliary file TFILE 
is require! , where entire Strip Tree corresponding to the 
feature on which computations are to performed, is copied 
from file STREE. This is essential since the records 
are required hot in sequential order of their existence 
in. the file, but as per the Tree structure, where pro- 
cessing is done in. order of LSON and RSON, recursively. 

The records are accessed more than once and appropriate 
ones selected in each sequence of processing. 
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4.6.3 Some Important Algorithms 

The basic operations on the Strip Trees were dis- 
cussed in Chapter XIX. Essentially all the computation 
to answer various queries make use of one or more of these 
basic operation. Certain queries which do not involve any 
feature , but only the point or points, do not make use 
of Strip Trees. These are: 

(a) Distance of a point to another point 

(b) Mag, True or Grid bearing of a point to another 
point. 

The algorithns for various queries are briefly discussed. 

4 .6.3.1 Distance 

For computing distance of a point to a feature, at 
resolution w r , LIST of nodes is made as discussed earlier. 
Distance from point to each of the strips in the LIST is 
computed, and minimum of these distances is reported as 
the result. 

To compute the distance of- point to a strip, it 
is computed as perpendicular distance to the line B-E, 
if the point is in region 'A* of the strip (point PI) 
or minimum of the distance from point B and E, in case 
of point lying outside the region 'A* (point P2). 

(Refer to Fig. 4.2). 
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4 .6,3 *2 Length 

As discussed earlier in Section 3.4.2, the error in 
length calculation can be quite unpredictable if strips 
are used. To avoid this problem, the actual length of 
the curve, for the portion enclosed in its corresponding 
strip, is also stored in one of the fields of the strip. 

In order to calculate length, it is sunmed up from 
the LIST of strips. Only in case, where the length of 
portion of curve is required, a modification is required. 
The two points which describe the portion of feature will 
specify the two nearest strips from the LIST. The length 
from these two strips to the next strips on both sides is 
computed as the shortest distance between the point and 
the next strip (Fig, 4.3), and adding the length of remai- 
ning strip as explained earlier. The length of curve as 
shown in Fig.- 4.3 will be - 

length « d.^ + length information of strips s 2 + 

length information of strip S 3 + ... + d 2 

where is minimum distance from point P x to 

strip S 2 

do is minimum distance from point P 2 to 
2 

strip S 4 . 
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4. 6. 3. 3 Intersection ( Crossing ) 

It requires tree structure to compute whether two 
features intersect each others or not, Algorithm discussed 
in Section 3.4.3 is used. 

In case points of Intersection are also required, it 
is obtained by finding point of intersection between each 
strip of LISTl (corresponding to first feature) and every 
strip of LIST2 (corresponding to second feature). 

4. 6. 3. 4 Test for Containment 

Algorithm of Section 3.4.4 is used. Semi-infinite 
line S Q is formed by chosing other point as (100.0, 100.0). 
Number of intersections is calculated as explained in 
Section 4, 6. 3.3, and result obtained, depending upon whe- 
ther this number is odd or even. 

4. 6. 3. 5 Bearing 

'Grid Bearing* from Point Pi to Point P2 is defined 
as, the inclination of vector CPT ~ p£ ) with vertical grid 
lines on the map. Thus 

Bearing (P 2, PI) * Bearing (Pi, P2) + 180° 

variations of Magnetic Bearing and True Bearing is. speci- 
fied on each map. Hence these can be computed from Grid 

Bearing. 
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4. 6. 3. 6 Height and Intervisibility 

Height of a point can be computed by testing its 
containment within the set of contours on the map. The 
height will be given by the contour with maximum height, 
which contains this point. 

The method involves testing containment of point, 
and once containment is found, to extract the information 
about height, of the contour. Presently, the implementa- 
tion of contours in creation of map is not organised 
to provide the height information and therefore the compu- 
tations for height can not be performed, once such an 
information is associated with contours, it can be imple- 
mented* However processing on HEIGHT command, and testing 
of argument has been implemented, and control is passed to 
appropriate Buffer procedure, which simply gives a message. 

Similarly computation of intervisibility involves 
computation of heights of the two points. In addition, 
height of highest point on the line joining two points is 
required • This has also been implemented only upto Buffer 
procedure level. 

4. 6. 3. 7 Area of a Closed Feature 

Like length, area of curve is also computed at the 
time of creating Strip Trees in module CREATTREE. This 
Is Stored in the identification line of each feature (In 
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case of open features the values is 0.0), To answer 
this query in 'Query Process' module, this value is obtai- 
ned from the identification line of the feature. 
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CHAPTER V 

CASE STUDY AND CONCLUSION 


5.1 Case Study 

The package has been tested using the base map £2^# 
for all the queries except HEIGHT and INTERVISIBILITY, 
for the reasons explained in Chapter IV. Fig. 5.1 shows 
the map used for this purpose. Some examples about the 
nature of queries and corresponding answers# are shown 
in Appendix 'F*. Fig. 5.2 shows the Strip Tree Struc- 
ture for RIVER-1 at different resolutions. 

5.2 Summary 

The present work is a general purpose package which 
can be used for geometrical query processing where 
curves are represented by set of points# instead of 
mathematical relations. The approximation using strip 
Tree representation# makes the computation process effi- 
cient where a large set of points can be approximated 
to a small member. Some of the queries e.g. containment, 
intersection, where the answer is negative Ci.e. point 
is not contained in a curve or the curves do not interact) 
is resolved by first root node itself. 
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5*3 Suggestions 

(1) Presently creation of base map [ 2 ], answering 
of non-graphical queries [ 4 J , and this work 
has been done independently as separate pack- 
age. Integrating it into ..one- package will 
Increase its usefulness. 

(ii) As discussed earlier, the contours on the base 
map generated, do not have the height infor- 
mation associated with them. If the necessary 
modification is made to include height infor- 
mation, queries of HEIGHT and INTERVISIBILITY 
can also be resolved. 

(iii) At present, answers to the queries are given 

in alpha-numeric form only, and they can not \ 

be displayed graphically. This can be made I 

f 

possible with raster graphics devices, where 
selective display and eraser is possible. ; 

whenever the answer to a query is some feature, [ 

I 

it can be indicated on the displayed map using f 

1 

f 

flashing arrows. | 

(iv) The concept of Strip Trees can be extended to 1 

f 

three dimensions to approximate heights, using \ 

t 

Strip Boxes. This may be useful in 3-dimen- i 

sional Sand Model displays. \ 

\ 

i 

i 

\ 

1 

I 

[ 
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(v) The package has used sequential files# where 
a sequential search for data pertaining to 
the required feature is not the most efficient 
solution. Suitable memory management using 
lower level language# and organising random 
access data files# the package can be made 
more efficient* 
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APPENDIX *A* 

TYPE OF GRAPHICAL QUERIES 

1 . Distance from a point to another point or a feature 

2. Length of a given feature 

3 „ To test whether two features intersect each other 

4. To test whether a point is contained inside a closed 
feature 

5 . Area of a closed feature 

6 . Height of a given point from contour information on the 

map | 

7 . To test whether two points are intervisible to each 

i 

others . s 

8. To find Magnetic , Grid or True bearing between two | 

f 

points . j 
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APPENDIX ‘C* 

TYPES OF ERRORS DETECTED 

1 . Illegal query type . 

2. Illegal argument, 

3. No argument to query, 

4. No HELP available on specified query, 

5. Too few arguments for this query, 

6 . GR contains non-integer characters, 

7 . GR is not a 6-figure GR. 

8. Illegal feature name. 

9. Feature name incomplete. 

10, Given feature name does not exist on map under refe- 
rence . 

11. Both arguments are same point - Bearing calculation 
not possible* 

12* Argument for this query should be a closed feature. 

13 . Both arguments for this query should be points . 

14 . Arguments for this query should be - A point and a 
feature name • 

15. Arguments for this query should be - A point and a 

i 


closed feature 
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16. Only argument for this query should be - point. 

17. Both arguments for this query should be - feature 
names. 

18. Argument for this query should be - closed feature. 

19. For LENGTH - if two arguments are points third should 
be feature name. 

20. Ambiguous INSIDE command. 

21. Ambiguous ENTIRE command. 

22. Illegal arguments for LENGTH. 

23. Illegal arguments for BEARING. 

24. Computations not possible. 
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APPENDIX *D' 


BNF GRAMMER FOR COMMAND LANGUAGE 


COMMAND : := <( UUERYNAMe) <( aRg) <^ARg) <)\Rg) 

QUERYNaME ::= <^HELp) | <HELp) /( . ^QUERYTYPE) 
^ QUERYTYPE ) 


QUERYTYPE : ^DISt) \ LENGTH) \ <CROSs) | -<AREa) 

^ CONTAIN) ] ^BEARING ) | ^HEIGHT) | 
^INTERVISIBILITy) | ^DISPLAY) 


ARG ss* >j ^ARGUMENT) 

ARGUMENT : := ^POINt) | ^FEATURECLASs) <^. ) <(eXTn) J 
^ENTIRE ) \ ^INSIDE) | <MAg) | <^GRID ) j 

<(true ) 

POINT < DIGIT) < DIGIT) <DIGIt) <Wt) < DIGIT )■ 

<^DIGIT ) 

DIGIT <(0) \ <1> | <2 > J <3 > | <4 >l<5 >|<6> l 7 

< 8 > i< 9 > 


FEaTURECLASS 


t= ^ OPEN ) \ <(CLO SED ) 

OPEN : := <ROADNh) I <ROADSh) [ <CARTTRACK> | <RIVER> 
<^railline) | <intboundary)> I <canal) 


CLOSED 


^contour) | <(enocarea) j ^ownocarea) j 


^FOREST ) | <^4ARSh) 
EXTN ::= <GlGIT) <^DIg)<DXG^> 
DIG ::= < J\ ^DIGIT^ 



APPENDIX 'E* 


ORDERED LIST OF QUERIES, FEATURE CLASSES AND 
SPECIAL ARGUMENTS 


( a ) Queries 

Ser. No « 

1 . 

2 , 

3 . 

4. 

5. 

6 . 

7. 

8 . 

9. 


Query Name 

DISTANCE 

LENGTH 

CROSSING 

AREA 

CONTAINMENT 

INTERVISIBILITY 

HEIGHT 

BEARING 

DISPLAY 


(b) Feature-Classes 


Ser * No * 

Ctoen Features 

Closed Features 

1. 

ROADNH 

CONTOUR 

2. 

ROADSH 

ENOC-AREA 

3. 

CARTTRACK 

OWNOC-AREA 

4 * 

RIVER 

FOREST 

5. 

RAILLINE 

MARSH 

6 . 

INTBOUNDARY 


7. 

CANAL 




